Pairing symmetry conversion by spin-active interfaces in superconducting junctions 
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We study the proximity-induced superconducting correlations in a normal metal connected to a superconductor 
when the interface between them is spin-active and the normal metal is ballistic or diffusive. Remarkably, for any 
interface spin polarization there is a critical interface resistance, above which the conventional even-frequency 
proximity component vanishes completely at the chemical potential, while the odd-frequency component re- 
mains finite. We propose a way to unambiguously observe the odd-frequency component. 
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Superconductivity and superfluidity are hallmarks of the 
wave-like character of matter, and manifest themselves in 
vastly different systems, from ultracold dilute gases via cold 
metals and fluids, to extremely dense protonic and neutronic 
matter In all these contexts, the symmetry of the order param- 
eter is of profound importance. Over the last decades, the pos- 
sibility of superconducting order parameters that change sign 
under a time-coordinate exchange of the two fermions com- 
prising the Cooper-pair, has emerged in addition to the by now 
well studied varieties of orbital symmetries IL^eiSe AtS,]. This 
so-called odd-frequency superconductivity ||6|] is distinct from 
the traditional even-frequency pairing in the Bardeen-Cooper- 
Schrieffer paradigm, and may be induced by proximityeffects 
in hybrid structures of superconductors and magnets Ij]] . 

In a broader context, proximity systems offer the possibility 
of controlling the physics of competing broken symmetries. 
The fundamental heterostructure for studying proximity in- 
duced superconductivity is the superconductor/normal metal 
(S|N) bilayer, where the normal metal or the interface may 
have magnetic properties. Among possible triplet pair coiTe- 
lations, in the diffusive limit odd-frequency pairs are favored 
I?], whereas in ballistic hybrid systems both odd- and even- 
frequency amplitudes compete |l3l|4|]. As all known supercon- 
ductors to date exhibit an even-frequency order parameter, the 
observation of proximity induced effects that are particular to 
odd-frequency pairing would be of utmost interest. 

There are two major difficulties associated with the detec- 
tion of the odd-frequency state in superconductor/ferromagnet 
(S|F) bilayers. One is the usually short penetration depth 
into the ferromagnetic region, limited by the magnetic coher- 
ence length £,F, much less than the superconducting coherence 
length iH . Another problem is that odd-frequency pairs are 
only well defined when even-frequency correlations vanish in 
the ferromagnet. Clear-cut signatures of the former are there- 
fore only accessible in a limited parameter regime |8]. 

The majority of work on superconducting proximity- 
structures so far has been restricted to the diffusive limit and 
spin-inactive interfaces |9]. For a non-magnetic bilayer, a 
minigap appears in the density of states of the normal metal. 
It scales with the Thouless energy of the normal layer and 



with the transmission probability of the interface. Such mini- 
gap structures are readily accessible experimentally ifioll . For 
a spin-active interface, the transmission properties of spin-f 
and spin- J, electrons into a metal are different, and this gives 
rise to both spin-dependent conductivities and spin-dependent 
phase shifts at the interface El El HH El- In this Letter 
we show that a spin-active interface in a S |N bilayer produces 
clear signatures of purely odd-frequency triplet pairing ampli- 
tudes that can be tested experimentally. 

We consider the system shown in Fig. [T] The superconduc- 
tor is conventional (even-frequency s-wave) while the inter- 
face is magnetic. We find that there is a dramatic change in 
the nature of proximity correlations when the spin-dependent 
phase shifts exceed the tunneling probability of the interface. 
The spin-active interface in an S|N bilayer causes the even- 
frequency correlations to vanish at zero excitation energy, 
while odd-frequency correlations appear At the same time, 
the minigap, one of the hallmarks of the conventional prox- 
imity effect, is replaced by a low-energy band with enhanced 
density of states. We focus on the density of states (DOS) 
in the normal region, which can be probed by tunneling ex- 
periments. Our findings suggest that it should be possible to 
detect the odd-frequency amplitude without any interfering ef- 
fects of even-frequency correlations. Since the exchange field 
is absent in the normal metal, this resolves the two main dif- 
ficulties associated with the experimental detection of odd- 
frequency correlations mentioned above. 

We adopt the quasiclassical theory of superconductivity 
1 1811 . where information about the physical properties of the 
system is embedded in the Green's function. For equilibrium 
situations, it suffices to consider the retarded Green's func- 
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FIG. 1: (Color online) Proposed experimental setup for observa- 
tion of the odd-frequency component in a diffusive normal metal 
layer I superconductor junction. 




2 



tion, g, that is parameterized conveniently in the normal (N) 
region by a parameter 0^, allowing for both singlet and triplet 
correlations |8]. In the superconducting (S) region, we em- 
ploy the bulk solution 55 = c • ts ® ctq + s • ri 55 (1(72), with 
c = cosh(0), s = sinh(0), 9 — atanh(A/£), and ai being 
Pauli matrices in particle-hole and spin space, respectively. 

We use the formalism described in Ref. |8], and con- 
sider first the diffusive limit. Then, the orbital symme- 
try for all proximity amplitudes is reduced to s-wave and 
hence the singlet component always has an even-frequency 
symmetry while the triplet component has an odd-frequency 
symmetry. The Green's functions are subject to boundary 
conditions, which in the tunnelinglimit assume the follow- 
ing form at the S|N interface |13l [Tsl]: 2^dgNdxgN ~ 
[ffSj 5Ar] + i(G0/G't)[2o "X'^jSa']; and at the outer interface 
read dxfjN = 0. Here, 7 = Rb/Rn where Rb (Rn) is the 
resistance of the barrier (normal region), and d is the width of 
the normal region, while Gt is the junction conductance in the 
normal-state. The boundary condition above contains an addi- 
tional term compared to the usual non-magnetic boundary 
conditions in Ref. |9]. This term is due to spin-dependent 
phase shifts of quasiparticles being reflected at the interface. 
Gcf, may be non-zero even if the transmission Gt — ^ 0, corre- 
sponding to a ferromagnetic insulator 1 13]. We define the su- 
perconducting coherence length ^5 ~ y^D/A and Thouless 
energy exh = D/d^, where D is the diffusion constant, and 
assume that the inelastic scattering length, Im, is sufficiently 
large, such that d <C ^in- 

The Usadel equation d reads Ddl9a + 2i£sinh6'„ = 0, 
with boundary condition ^ddxOa = {cSa — crsca) + lasa^ 
at X = and dxda- — al x = d. Here, Ca — cosh(6'o-) 
and Scr = sinh(f?o-). At zero energy, we find that the pairing 
amplitudes are either purely (odd-frequency) triplet. 
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Thus, the presence of G^ induces an odd-frequency compo- 
nent in the normal layer. The remarkable aspect of Eqs. (HJ 
and (I2J is that they are valid for any value of the width d be- 
low the inelastic scattering length, and for any interface pa- 
rameter 7. Thus, the vanishing of the singlet component is 
a robust feature in S|N structures with spin-active interfaces, 
as long as jG^j/Gr > 1- Without loss of generality, we 
focus on positive values of G^ from now on. The DOS is 
given as N{e)/No = E.Re{c.}/2, yielding A^(0)/iVo = 
RejG^/^G^ - G^}. At zero-energy, the DOS thus van- 
ishes as long as G^/Gt < 1, which means that the usual 
minigap in S|N structures survives in this regime. However, 
the zero-energy DOS is enhanced for G^ / Gt > 1 since the 
singlet component vanishes there. 
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FIG. 2: (Color online) The singlet and triplet proximity amplitudes 
induced in the normal metal are shown for G4, / Gt < 1 [in a) and c)] 
and G^/Gt > 1 [in b) and d)]. In e), we plot the energy-resolved 
DOS for several values of G4, / Gt- Finally, f) shows the zero-energy 
DOS as a function of Gt / G^, with the proximity amplitudes shown 
in the inset. 



The full energy-dependence of the DOS may only be ob- 
tained numerically. To model a realistic experimental setup, 
we fix 7 = 10 and d/£^s = 1-0, although our qualitative re- 
sults are independent of these particular choices. As a mea- 
sure of the relevant energy scale, we define Eq = eTh/(27). 
The results are shown in Fig. |2] to investigate the effect 
of the spin-dependent phase shifts. The low-energy DOS is 
strongly enhanced due to the odd-frequency amplitude when 
G^/Gt > 1 (G^/Gt = 1-5 in the figure). Conversely, the 
DOS develops a minigap around e = when G^/Gt < 1 
(G^/Gt ~ 0.5 in the figure). The ratio G^/Gt depends on 
the microscopic barrier properties l"!^. In the tunneling limit, 
one finds that G^ can be considerably larger than Gt- 

We suggest the following qualitative explanation for the 
mechanism behind the separation between even- and odd- 
frequency correlations. The superconductor induces a mini- 
gap oc Gt in the normal metal, while the spin-active barrier 
induces an effective exchange field oc G^. The situation in the 
normal metal then resembles that of a thin-film conventional 
superconductor in the presence of an in-plane external mag- 
netic field i22ll . with the role of the gap and field played by 
Gt and G^, respectively. In that case, it is known that super- 
conductivity is destroyed above the Clogston-Chandrasekhar 
limit 1 23], as the spin-singlet Cooper-pairs break up. In the 
present case, we observe coexistence of the exchange field 
and spin-singlet even-frequency superconductivity as long as 
Gcf, is below the critical value of G^ — Gt- However, for 
Gcf, > Gt spin-singlet pairing is no longer possible at the 
chemical potential. It is then replaced by spin-triplet pairing, 
which must be odd in frequency due to the isotropization of 
the gap in the diffusive limit. Thus, there is a natural sepa- 
ration between even-frequency and odd-frequency pairing in 
the normal metal at a critical value of the effective exchange 
field G^. 

The same effect occurs in the ballistic limit, as we now 
show. In this case, we can obtain the retarded Green's func- 
tion using the formalism described in Refs. Hill. The 
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Eilenberger equation in the normal region reads ivp^dxg + 
[£T3 ® (To, g] =0. For the boundary conditions, we use a scat- 
tering matrix describing the magnetic interface between the 
superconductor and the normal metal, 



S 



Ins ■ exp {^'&ns^ 



isN ■ exp (fl^SAfOs) 

-Ln ■ exp i^i^NOa) 



(3) 



with real reflection and transmission spin matrices r^, r^, 
tg]^, and tpfg. The spin mixing angles iJs, ^n, "dsN, and 
-djsis describe spin dependent scattering phases Ml ill . Neglect- 



ing spin flip scattering, the transmission and reflection ampli- 
tudes are diagonal in spin space, and the relations rg = Zn = 



diag [rT,rx],i^5 = 



diag [t|,tj,r^-ft2 



1, 



"^NS + "^SN = "^s + '^N follow from the unitarity of S. Pos- 
sible scalar phases are omitted in Eq. ([3]l, as they play no role 
in the final results. 

We next concentrate on subgap energies. The anoma- 
lous amplitudes can be decomposed into singlet and triplet 
components, / = + fto^jiiaj). Defining fa = 
(/s + a ft)/ 2, we obtain on the top of the normal over- 
layer (a; = d) fa{e) = -sgn{aa)t^ti/y/al - {t^tiy wifli 
Qfo- = siTi{2sd/vpx + a"S+) + r^rism{2£d/vFx + o''d-). 
Here, 'd± = ^{-diy ± -ds) ± arcsin(e/A), and e has to be 
supplemented by an infinitesimally small positive imaginary 
part. The interface parameters and the Fermi velocity com- 
ponent in x-direction, vpx — vp cos4>, depend on the impact 
angle tp. The relevant energy scale in the problem is the bal- 
listic Thouless energy, eth — vp/2d. For zero spin mixing 
angles we recover the known DOS for a normal state over- 
layer on a singlet superconductor The DOS is non-zero only 
for |a| > which for sufficiently large impact angle al- 
ways is fulfilled. Clearly, the most interesting regime concerns 
e/sTh |i?±| t^ti. 

In the tunneling limit, for small excitation energies 
s/sTh ■C 1 and small spin mixing angles d± we obtain 
Q!o- = {4:ed/vpx + (J'&n)- In this case, due to + — i?jv, 
only the spin mixing angle for reflection at the normal side 
of the interface enters, and acts as an effective exchange field 
b — "df^ypx/idonthe quasiparticles. Especially interesting is 
the case £ = 0, for which all proximity amplitudes are even in 
momentum. For e = we obtain aa — cri9jv, and the pairing 
amplitudes are either purely (odd-frequency) triplet. 
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Comparing with the results for the diffusive case, we find that 
G^/Gt corresponds to — i?Ar/(tf ij). 

In Fig. |3] we show results for the proximity amplitudes in 
the ballistic limit, and focus on positive values of-dM without 
loss of generality. A systematic expansion of all terms in the 
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FIG. 3: (Color online) Momentum-averaged proximity amplitudes at 
the surface of the normal layer. Parameters: d — vp/A, To = 0.1 
(see text), (a) and (c): j?iv = j?s = 0.05 < To; (b) and (d): 
i?jv = = 0.15 > To. Energy units are eo = To erh- Even fre- 
quency singlet components are shown in (a-b), odd frequency triplet 
components in (c-d). 












9^=0.05 





9^^=0.15 



1.5 
1 

0.5 


-0.5 



(b) j^(Q)/% 8„=o.2 




0.2 0.4 0.6 0.8 1 







e 


=0 


— 6 


=0.05 


-- e, 


=0.15 



^0. 



(d) 



0.02 
0.01 


-0.01 
-0.02 



ImffOJ _ 



Re/fOJ z 



0.5 



0.2 0.4 0.6 0.8 1 







FIG. 4: (Color online) (a) DOS as function of energy at the top of 
the normal layer for fixed transmission probability To — 0.1, and 
various values of i^at = i?s. Remaining parameters are as in Fig. [5] 
(b) DOS and proximity amplitudes at e = for '!9jv = i9s=0.2as 
function of To. In (c) and (d) we show the results corresponding to 
(a) and (b) when assuming an (abrupt) tunneling cone with opening 
angle of 10 degree. 



tunneling probability shows that in the tunneling limit the spin 
dependence of the transmission probabilities can be neglected, 
and only that of the phase shifts needs to be kept. Thus, we 
assume ~ ~ t. We model the dependence on the impact 
angle xjj as = (<o) ^ M = cos-0, and assume for simplic- 
ity spin mixing angles independent of ii. The tunneling prob- 
ability for normal impact is Tq — t^. In the case To < iSn at 
small energies the odd frequency triplet amplitude dominates, 
and it is the only non-zero amplitude at e = 0. On the other 
hand, for Tq > -d^ both singlet and triplet amplitudes con- 
tribute. This is due to the fact that for large impact angles the 
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transmission probability drops below the value for the 
spin mixing angle i^n ■ 

We turn now to the DOS. The general expression, assum- 
ing the bulk solution in the superconductor, is N{e)/NQ = 

ReEa=±i lo \^'y\N^l - ^h^iY d^i. In the tunneling 
limit, this simplifies again, and provided that |i?Ar| > 
i|t| for all impact angles, the DOS at the Fermi level 
is enhanced above its normal state value, N{0)/No = 
/ dfi I'^^nI/ — (t|tx)-^. In Fig. m we show results for 
the DOS. In (a-b) we assume the dependence on the impact 
angle as above, whereas in (c-d) we allow tunnelling only in a 
narrow tunneling cone of 10 degrees. The DOS for the cases 
of dominating triplet amplitudes and dominating singlet am- 
plitudes differ qualitatively. In the case of a tunneling cone 
this difference is most drastic, and a comparison with the re- 
sults above shows that it is very similar to the diffusive case. 
In the right panels, where — "ffs — 0.2, we demonstrate 
that for Tq < 0.2 only the odd frequency triplet amplitude 
is present at the chemical potential, while the singlet ampli- 
tude is zero. The corresponding zero-bias DOS is enhanced 
in this region, whereas it is reduced in the region when singlet 
correlations are present at e = 0. 

The simplest experimental manifestation of the odd- 
freq uency component is a zero-energy peak in the DOS 
1 17, l20i u\\\ . In S|F layers, where this phenomenon has 
been discussed previously, a clear zero-energy peak is unfor- 
tunately often masked by the simultaneous presence of singlet 
correlations fg, which tend to suppress the DOS at low ener- 
gies. This is not so in the system we consider, provided only 
To < |i?Ar I in the ballistic limit, or equivalently, Gt < jG^jin 
the diffusive limit. This is ideal for a direct observation of the 
odd-frequency component, manifested as a zero-energy peak 
in the DOS. 

The important factor, with regard to isolation of the odd- 
frequency correlations at zero energy is the interface. The 
even-frequency correlations vanish completely when the in- 
terface transmission is sufficiently low. The parameters ^n, 
or equivalently, G^ can be increased by increasing the mag- 
netic polarization of the barrier separating the superconduct- 
ing and normal layers. By fabricating several samples with 
progressively increasing strength of magnetic moment of the 
barrier, one should be able to observe an abrupt crossover at 
the zero-energy DOS above a certain strength of the magnetic 
moment. Alternatively, one could alter the interface transmis- 
sion by varying the thickness of the insulating region. 

In summary, we have investigated the proximity-effect in 



a superconductor/normal metal bilayer with spin-active inter- 
face. We find that both in the ballistic and diffusive limits, 
the usual even-frequency correlations may vanish completely 
at zero energy, while odd-frequency correlations persist. This 
result is completely independent of the specific values for the 
layer thicknesses and barrier resistances, indicating that it is a 
robust and general feature of spin-active interfaces. Our find- 
ings suggest a way of obtaining unambiguous experimental 
identification of superconducting odd-frequency correlations. 
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